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The use of operator methods of algebraic nature is shown to be a very powerful tool to deal with 
different forms of relativistic wave equations. The methods provide either exact or approximate so¬ 
lutions for various forms of differential equations, such as relativistic Schrodinger, Klein-Gordon and 
Dirac. We discuss the free particle hypotheses and those relevant to particles subject to non-trivial 
potentials. In the latter case we will show how the proposed method leads to easily implementable 
numerical algorithms. 


I. INTRODUCTION 

Operation techniques embedded with the formalism of 
the evolution operator have played a central role in the 
study of time dependent problems in quantum mechan¬ 
ics. The associated technicalities have attracted the at¬ 
tention of physicists [I[ and mathematicians 0 , who pro¬ 
vided tools of crucial importance to deal with the study 
of the time dependent Schrodinger equation. Expansions 
of the evolution operator like those developed by Stueck- 
elberg, Feynman and Dyson on the physicists’ side and 
by Magnus and Fer on the other have paved the way to 
the formalism underlying the modern diagrammatic pro¬ 
cedures il. Evolution problems in relativistic quan¬ 
tum mechanics have raised new questions mainly associ¬ 
ated with the pseudo-differential nature of the involved 
operators. These methods, neither popular nor well ap¬ 
preciated by physicists in the past, are now getting in¬ 
creasing attention for a more sound formulation of the 
underlying mathematical foundations and for the wealth 
of phenomenology they can describe In particular, 
the operational approach has been shown to be effective 
in solving pseudo-differential equations, as well as those 
containing fractional derivatives and related to anoma¬ 
lous transport phenomena. The latter include not only 
the so-called fractional dynamics adopted to investigate 
anomalous diffusion and similar effects in complex phys¬ 
ical systems Q but appear useful to model problems in 
life sciences, engineering, economy, and even studies of 
human mobility [Ml. 

In this paper we will focus our attention on the equa¬ 
tions of the relativistic quantum mechanics. Our aim is 
to push forward previous investigations of solutions to the 
relativistic evolution equation involving fractional deriva- 
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fives, namely the relativistic Schrodinger (called also 
Salpeter’s) equation [H, [11, [il|, and to compare them 
with operationally obtained solutions to the Dirac-type 
and Klein-Gordon equations. The main motivation of 
this research is that fundamental relativistic wave equa¬ 
tions give results whose interpretation based on Hamilto¬ 
nian mechanics concepts, as well as explanations of aris¬ 
ing paradoxes (like Zitterbewegung or the Klein para¬ 
dox), still remains controversial [l7| and need clarifica¬ 
tion. In our opinion the use of fractional derivatives or of 
pseudo-differential operators, to treat relativistic quan¬ 
tum problems, allows a more complete understanding of 
longstanding problems associated with nonlocal evolu¬ 
tion problems. These problems, once considered aca¬ 
demic, have in recent years acquired a more concrete 
flavour, mainly in the context of the so-called quan¬ 
tum simulation, namely the possibility of ’’simulating” 
quantum nonrelativistic effects in condensed matter and 
atomic physics (l^ . 

We begin with recalling results coming from the study 
of the relativistic heat equation (as proposed in O) 
whose solutions provide us an example of the method. 
The Section III is devoted to the analysis of the free rela¬ 
tivistic Schrodinger equation and unexpected behaviour 
of its solutions. In the Section IV we consider the Dirac 
factorization of the relativistic Schrodinger equation and 
compare so obtained results with those of Section III. 
The Klein-Gordon equation, together with the relation 
of its solutions to the Dirac case, is studied in the Sec¬ 
tion V, while in the Section VI we present the analysis 
of the simplest “interacting model”, namely the relativis¬ 
tic particle (Schrodinger’s and Dirac’s) driven by a linear 
potential. In what follows we restrict ourselves to 1-1-1 
spacetime models. It simplifies the presentation and pre¬ 
serves readability of the results but there is no reason to 
assume that general properties of solutions will not be 
shared with the realistic 1-1-3 situation. 
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II. THE RELATIVISTIC HEAT EQUATION • The integral kernel 


In Ref. [l^ the following equation has been defined as 
the relativistic heat equation 

dtF{x,t) = —^Jd? — k'^d^F{x,t), (1) 

F{x,0) = f{x), (2) 

where t has the dimension of time, d and k have the 
dimensions of the inverse time and of a length divided by 
time, respectively. Introducing dimensionless variables 
T = dt and ^ = xd/k we can cast Eq. ([T]) in the form 


drF{i,T) = -^l-dlF[^,T). (3) 

From the mathematical point of view Eqs. (Ill) and ® 
are evolution equations with the distinctive feature that 
they contain the square root of a differential operator. 
The formal solution of Eq. can be written as 

F{^,T)=e--V^f{0- ( 4 ) 


Albeit the use of the exponential operators is quite 
common, the handling of exponentials containing the 
square root of a differential operator (as well as higher 
order roots) have been rarely treated. Nevertheless, some 
analytical tools have been developed for this purpose 
in the past. In the specihc case of Eq. O a suitable 
representation, which allows evaluation of the action of 
’’square rooted” operator on the function /(^), can be 
realized through the Doetsch formula . The latter al¬ 
lows explicit evaluation of the r.h.s. of Eq.Q according 
to the following integral transform 





9i/2iv)e 


fiO, 



5i/2(??)e ^^^D{^,r]T^)dr]. 


( 5 ) 

( 6 ) 


where D{^,r]T^) = e’'” « fiO- Note two peculiar fea¬ 
tures of the Eq. ®: 

• Solution of the relativistic heat equation Eq. ® 
depends, via an appropriate transform, on the so¬ 
lution of the “classical” diffusion problem as D{^, r) 
is the solution of the ordinary heat equation 


drD{^,T)=dlD{^,T), (7) 

D{^,0) = f{0, ( 8 ) 


usually expressible in terms of the Gauss- 
Weierstrass transform 

OO 

D{C,t) =e^^(f{C) =-^= [ e-^^^f{a)da. (9) 

2V7rr J 

— OO 

or written, using the Fourier transform, as 

1 

= ^ F‘>^e^‘>^f{q)dq, (10) 

V J —CO 

with f{q) = (2^)-i/2/-^e-*««/(e)de 


is a generic example of the so-called Levy stable 
distributions or a-stable distributions. The latter 
denoted as ga(,q) are for r; > 0, and for 0 < a < 1, 
given uniquely by the inverse Laplace transform of 
i.e. e~P°‘ = e~P^gaiq)dr] see [2ll-l^. 

Let us now consider a solution of the relativistic heat 
equation for two specific initial functions. First, we as¬ 
sume that the initial condition of Eq. ([5]) is a Gaussian: 


fmo = 


2^/Tr^ 


exp 


4/3, 


h 




( 12 ) 

where /3 > 0. We can therefore write the relevant solution 
as 


-Di(/3;Ct) 


1 

— , exp 


e 

4(/3 + t) 


1 

J — CO 

recognized as the so-called Glaisher identity 
Eq. (fT^ . if substituted into Eq. ®, leads to 


(13) 



Fi{I3-,^,t) 


1 poo 

27 r J_oo 


(14) 


the same as the result of formal manipulation on Eqs.® 
and (fl^ . 

The second example of the initial condition under con¬ 
sideration is 

1 fOO 

= mm I-F 


where /3 > 0, Ki{j3) is the modified Bessel function of 
the second kind, i.e. the Macdonald function [2^. Eq. 
(fT51l is essentially the Fourier transform of 


/2(/3;fc) 


vT+F- 

2Ai(/3) 


(16) 


normalized in the £^-norm: f2{/3; k)dk = 1. The 

solution of the relativistic heat equation with Eq. m 
as initial function, calculated as shown before, reads 


F2{P;^,t) = 




/ 

J —( 




277i(/3) J-c 

1 poo _ 

The use of the identity (see Refs. [H, [2^) 

[ Qos(&x)(ia: = 2aq - K^^qy/o? 

7-00 ’ 

(18) 
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allows to cast Eq. (II3 in the form 


P + T Ki{^{I3 + tY 


2^1 (/ 5 ) ^{P + rY+e 


which may be viewed as the identity 


^ P + T Kl{y/{P + TY +C^) 

" 2 Ki(^) ^(/3 + r) 2+^2 


(19) 


( 20 ) 


and interpreted as a generalization of the Glaisher for¬ 
mula for the Macdonald function. 

Concluding this section we would like to note that 
naming Eq. m as relativistic heat equation may lead to 
some misconception and should be used carefully, or even 
in quotes. The concept of relativistic diffusion needs clar¬ 
ification at least from the mathematical point of view. 
Namely Eq. o does not represent a diffusion process in 
the strict sense because the norm of the relevant solu¬ 
tion is not preserved during the evolution and decays 
with 6“’’. Eqs. da and dl) are properly normalized 
at any time if multiplied by and Ki{P)e^^'^/n, re¬ 
spectively. In such a case they may be considered as the 
evolution of a distribution which satisfies the equation 
3rEi(^,T) = -[(1 - a|)^ - l]Fi(^,T), i = 1,2 and ex¬ 
hibits the correct nonrelativistic limit in the sense that 

-[(l-a|)i-l]c.ia|-|9| + .... 


III. THE RELATIVISTIC SCHRODINGER 
EQUATION (SALPETER EQUATION) 

The example of relativistic evolution, which will be 
considered here, is the free relativistic Schrodinger equa¬ 
tion, called also the Salpeter equation, which has the 
form 

ihdt'^ = c^Jp'^ + moC^di'. (21) 

Written in dimensionless variables it reads 

idr'^{q,T) = y/l+^^{q,T), (22) 

where r = ctjXc, q = Xcp/h, and Ac = h/{mc), and then 

*a,vI/(e,r) = ^l-5|vI/(e,T), 

= (23) 

with q = id^, ^ = x/Xc, originating from the standard 

quantization rule p = ihdx- 

In analogy to the Eq. ([51) we can express the formal 
solution of Eq. (051) in terms of a Gaussian wave function 
D{^,t) and the Levy distribution gi/ 2 {'q) 

vl/(^,r)=e-*"v^vI/o(^) 

poo 

= / 9i/2{v)e'^'"^ D{^,-r]T'^)dT]. (24) 

Jo 


Taking D{^,—r]T'^) in the Fourier space (see Eq. (ITOl) 
with appropriately changed variables) we can show that 
Eq. (l24ll is equivalent to the solution of Eq. (l2^ . It is 
seen from 


poo r poo 

0 —oo 


dq 


drj 


I 




— CO 0 

1 


J-c 




^o{q) 


dq 

(25) 


The inverse Fourier transform of Eq. (|25|) defines ^{q, r) 
which is a formal solution of Eq. (051) . 

Let us now solve the relativistic Schrodinger equation 
for an initial Gaussian packet given in Eq. (El), i.e. for 
the initial condition e~^'^ taken in Eq. (1551) . Thus we 
have 

T\T poo _ 

vl/i(/3;^,r) = ^ / (26) 

V -J —oo 


with the normalization constant N = [2/?/(7rAc)]^/^ cal¬ 
culated in C'^\pZ,dx\ space. 

A time-dependent solution of Eq. (l23l) obtained along 
the same lines as before for To(g) = 2x7/3) 
reconstructs the result presented in [I^, see Eq. (4.18) 
there. Namely, we get 


^2(/3;C,t) 


7l(vZ±l!) 

V^Ki{2P) yW+F 
P + iT Ki{^y{P + irY +C^) 
^JnKipjp) ^{P + irY +SY 


A comparison between nonrelativistic and relativistic so¬ 
lutions reveals interesting behaviour illustrated in Fig. 
[T^) andllJo). The snapshots, taken at different times, 
show that the wave packet spreading is “slower” for the 
relativistic case than for its nonrelativistic counterpart. 
Furthermore, the relativistic solutions reveal, at longer 
times, the wave function deformation determining their 
departure from the Gaussian-like initial form. Numerical 
calculations of '^i{P] t) done for various values of P and 
T show that there exists a non-trivial relation between P 
and T which implies an appearance, or non-appearence, 
of two peaks. Two peaks evolution pattern appears for 
reasonable values of r if /3 is small enough; if not, the 
evolution follows the standard spreading. As a quantity 
indicating what is going on for a given r the value of 

the second derivative (/3 ,t) = d‘l\^i{P]iiT)p , 

i = 1,2, may be used. For r = 0 it is evident that 


Df'\p,Q) < 0 for both Eqs. (l26l) and (l27l) . But 
if T grows it may happen that D^^\p,t) changes its 
sign, which means that the point f = 0 becomes local 
minimum. This depends on the parameter /?, measur¬ 
ing the width of the initial wave packet, i.e. the ini¬ 
tial localization, as well. From numerical calculation of 
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initial condition defined in the Fourier space as follows 




jV|g| e-f^yA+^ 

(yrr^-1)1/2 ■ 


(28) 


After a simple calculation, applying Eq. (2.5.39.6) on p. 
456 of [ 2 ^, Eq. (1^ . with the initial condition Eq. (1^ . 
reads 


4/3(/ 3 ;^,r)=iV(v^^^TF+ «)'/" ^7^5, ( 29 ) 

where a = f] + ir. The normalization constant N cal¬ 
culated in is given as N = y/n/Xc[Ko{2/3) + ^ — 
^/3Xo(2^)L_i(2^)-7r/3iFi(2^)Lo(2/3)]-i/2, where L^w) 
is the modified Struve function. Now we can find the ex¬ 
act relation between /3 and r for which two peaks appear. 
Namely, 


1-9^ N e~^P 

3 20F(/32 +t2)5/2 

X [r2(3-4/3)-/32(3 + 4/3)], (30) 

from which it is immediately seen that 

• D^^\p,0) < 0 for r = 0, 

• if /3 > 3/4 then D^^\j3,T) < 0 for all r, 

• if /3 < 3/4 then D^\p, r) becomes positive for r > 
r, = (/32(3 + 4/3))/(3-4/3), 

The behaviour of |4'3(/3;r)p for /? = 1/4, for which 
Te = 72/4, are illustrated in Eig. [2] in which the case 
T < Tc is marked by blue line, r = Tc is denoted by red 
line, and r > Tc is green line. 


FIG. 1. (a): Comparison, for /3 = 0.2 and r = 0.2, of the 
relativistic solutions llki)/!;v)!, i = 1,2 (Eqs. (27) and (26) 
for lines I and II, respectively) and nonrelativistic solutions 
F'i(/3;5,r), i = 1,2 (Eq. (19) and (14) for lines III and IV, 
respectively), (b): the same as in Fig. 1(a) but for j3 = 0.2 
and r = 2. 


= ■ 

T we can find 


5|l'I^.(/3;7r)|^ 


«=o 


~ 0, i = 1,2, for fixed 


that there exist Pc such that for P > Pc 


we can observe two peaks. For instance, D^\p, 7) ~ 0 

for Pc = 0.689 and Df\pj) ~ 0 for Pc = 1.411. All 
this allows us to make a conjecture: two-peaks evolu¬ 
tion pattern is possible only if the initial wave packet is 
localized sufficiently strongly, P < Pc and the “critical” 
localization Pc is related to the Compton wavelenght of 
the particle. 


To justify the previous assumption we will consider 
another solution of Eq. (1^51) . which is obtained for the 
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FIG. 2. The three various behaviours of |4'3(/3;C,'r)|^ for 
P = 1/4 and r = 0.2 (blue line), t = Tc = s/^jA: (red line), 
and T = 1 (green line). 
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Let us now come back to the consideration in dimen¬ 
sional coordinates x and t. According to the choice of the 
normalization of variables the all initial wave function is 
localized within the Compton wavelength of the particle 
itself. For example, Eq. (l27l) expressed in dimensional 
form reads 

PXc + ict Ki{^y/{l3Xc -I- icty + x'^) 

'i' 2 (x,t) = — , =- , -. 

y'XcTrKi{2P) ^ (/3Ac -f ictY + x'^ 

(31) 

The role of the parameter /3 appearing in Eq. dSB) is not 
trivial. It is an arbitrary constant associated with the 
intrinsic nonlocality of the problem under study and it 
defines the scale of delocalization of the particle within 
its Compton wavelength. 

The wave packets whose evolution has been reported 
in Figs. [T]are initially localized within a Compton wave¬ 
length. In the case of electronic neutrino assuming a mass 
of 2 eV we are talking about an uncertainty in position of 
the order of few microns. The case of massless particles 
could be obtained from Eq. (IHTl) by going to the limit 
Ac —> oo. However keeping such a limit without any fur¬ 
ther caution does not provide physical results. Eq. dSlD 
yields indeed a completely delocalized function unless we 
assume that in such a limit /3Ac remains a constant. Since 
(3 is an arbitrary scale factor, we can define b = jSXc- 
By such a redefinition the dependence on the Compton 
wave-length is a physical dependence and no more a scale 
factor. We can keep the limit safely and using the well- 
known asymptotic property Ki{z)\z^-^ ~ 1/z, we obtain 
the same result as quoted in Ref. 




12b 


b + ict 


TT {b + icty + x'^ 

It is now worth to recast Eq. (1321) in the form 


27r \_b + i{x + ct) b — i{x — ct)_ 


(32) 


(33) 


and note that it can be viewed as the solution of the 
equation 


= (34) 


This aspect of the problem may acquire a deeper physical 
consistence when studied in the presence of a potential, 
as will be discussed in a forthcoming paper. 


IV. THE DIRAC FACTORIZATION AND THE 
ONE DIMENSIONAL QUANTUM 
RELATIVISTIC EQUATION 

This section is devoted to the Dirac factorization of the 
one dimensional relativistic Schrodinger equation. This 
problem has recently raised a significant amount of in¬ 
terest because it may provide an important toy model to 
benchmark genuine quantum relativistic effect like the 
Zitterbewegung, in the so-called quantum simulation of 
Dirac equation [28l| . The equation can be written in the 
two component form [2^ as follows 

ihdt'^ = cap^ + (36) 

In the dimensionless coordinates r and ^ it can be ex¬ 
pressed by 


idr^ = —-I- (37) 

where a and /3 are Pauli-like matrices satisfying the iden¬ 
tities 


= 8 ^ = 1 = 



a8 + 8a = 0 (38) 


and numerically given as 


a = 



and 8 



We have denoted by J;. = j the two component vec¬ 
tor, replacing the wave function di*. The vector should 
not be confused with a spinor and Eq. (1^ is not the rel¬ 
ativistic counterpart of the Pauli equation since we are 
not dealing with a problem including the spin degrees of 
freedom. The reason for using two dimensional matri¬ 
ces is that we want Eq. (157)) to describe negative and 
positive energy states. 

The solution of Eq. dSZl) can be formally obtained, 
using the evolution operator method, as 


which is the massless limit of the relativistic Schrodinger 
equation. Eq. (IMl) does not depend explicitly on the 
Planck constant; such a dependence is however implic¬ 
itly contained in the arbitrary constant b. We find also 
that the r.m.s. value of the particle position is indeed 
provided by yix{8^ = \/b'^ + {ctY which is consistent 
with the constraint \/(p(t)^) > ^Ii8‘\jb’^ + (ct)^) (for ad¬ 
ditional comments see Ref. M)- The emergence of two 
peaks, when the particle is localized according to the pre¬ 
vious discussion, is a distinctive feature of the nonlocal 
nature inherent to the relativistic Schrodinger equation. 


J^(C,r) = U(T)3L(f,0) (39) 

with the evolution operator written in the form of the 
2x2 matrix 


U (r) = exp 


df —i 


AA(t)-I-.B( r) ^(t)94 \ 

\ A(t)5j -iA^r) + B{t) J ’ 


( 40 ) 
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where A{t) = sin(\/Ar)/\/A, B{t) = cos(\/At), and 
A = 1 —9|. This evolution operator is unitary and there¬ 
fore the Dirac wave function remains normalized at any 
time, namely 

/ CO poo 

-oo —oo 

(41) 

The use of the Fourier transform method allows one to 
obtain the solution of the free particle Dirac equation: 
from Eq. (1391) with U{t) given in Eq. (l40ll we get the 
evolution of the two components in the form 

'4’+{.£.,t) = j I^ ['0+ (fc: 0) + (fc: 0)] 

-|-cos(\/l + t) Vj+(A:,0)| (42) 

= J I ^ (fc: 0) - (fc: 0)] 

-I-cos(\/r-M)^r) ■!/)_(/(:, 0) [ (43) 

' V Stt 

where il)±{k^ 0) are the Fourier transforms of 'ip±{^, 0). In 
Fig. [3] we present the snapshots of the time evolution 
for |4 'i(^,t)P, calculated using Eqs. (|4T]) . (1421) and (l43l) . 
with the initial wave packet being a Gaussian containing 
both positive and negative energy components 

ii(fc,0) = e-^^-^ Q (44) 





FIG. 3. Evolution of an initial wave packet with both energy 
components as in Eq. dMj for d = 1 and different r values: 
T = 4.5 (blue line), r = 6.5 (red line), and r = 8.5 (black 
line). 

The presence of the negative energy term yields the be¬ 
haviour with some features not exhibited by the Salpeter 
equation. These effects are due to the superposition be¬ 
tween the two states, which induces also the characteris¬ 
tic trembling motion of the packet average position and 
yielding the “Zitterbewegung” [3l| , namely an oscillating 


motion around a linear trajectory. We remark that using 
Eq. (H31) with the initial condition possessing only the 
positive energy term which is given by Gaussian func¬ 
tion, we get the evolution which is not dissimilar from 
the one described with the relativistic Schrodinger equa¬ 
tion and from the behaviour obtained in Ref. 1^, using 
different solution techniques. 

Following [23 we can also find an exact solution for the 
one dimensional massive Dirac equation with the initial 
condition 




(45) 


According to Eqs. (fT8|) . (l42l) and taking the integral rep¬ 


resentation of (I + A:^) 
we obtain 


2'|-l/2 


i.e. 


v/T+F’ 


= f^e-^^^ds, 


I 


+ 


I 


Ko{c*) + 
—Kq{c) + 


P — IT 


iFi(c*) 


• IT 


Ki{c) 


and 


'0-(C,t) = 






iFi(c) - -Mc^) 


(46) 


(47) 


where c = \/{P + *t)^ -I- and c* is its complex conju¬ 
gation. 


V. ONE DIMENSIONAL KLEIN-GORDON 
EQUATION 

The free particle Klein Gordon (FP-KG) equation [I[ 
has played a fundamental role in the description of rel¬ 
ativistic spin 0 particles. It is well known that accord¬ 
ing to Feshback Villars transformation [s^, FP-KG is 
transformed in two component Schrodinger like equation 
which exhibits a Hamiltonian with unconventional prop¬ 
erties. In this Section we treat this transformation from 
a different perspective. We will look for a formal solu¬ 
tion of the free particle KG equation by treating it as 
a D’Alembert equation and next we will show that the 
solution of the I-dimensional free particle KG can be ex¬ 
pressed in terms of those of the Dirac equation. 

Let FP-KG be written in the dimensionless coordinates 
^ and T introduced in Sec. Ill as follows 

92vl/(^,r) = -(l-5|)vl/(^,r) (48) 

with the initial conditions 

'f'(e, 0 ) = 4/0(0, [^^r 4 /(OT)],^o = 4/1(0. ( 49 ) 

Using the formal procedure presented in the previous Sec¬ 
tions for solving Eqs. (HSl) and (H31) we can cast our prob¬ 
lem in the form of an ordinary second order equation with 
constant coefficients, namely, 

924/(00 = - 024 /( 00 , 0 = (1-92)1/2, (50) 
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FIG. 4. Evolution of the FP-KG solution, Eq. (1551) . at differ¬ 
ent times: r = 0.5 (I; red line), t = 1 (II; blue line), r = 4 
(III; green line), and r = 9.5 (IV; gold line). 


which is resembling that of a harmonic oscillator. There¬ 
fore, we postulate that the solution of free particle KG 
can be written in the form 

'^kg{^:t) = + e-*"*C2(e), (51) 

where Ci ,2 are linked to the initial conditions of Eq. 
by the identities 

Ci+C2 = ^o(e), = (52) 

Thus, inserting the explicit formulas on Ci and C 2 into 
Eq. dSU and making simple symbolic manipulations, we 
get 

= i[i«'/j_s(^,T) +2^'ij-s(C,-T)] (53) 

= cos(rl))d'o(C) + sin(ri?)I)“^d'i(^). (54) 

The function (^, r) is the solution of the (forward) 

free relativistic Schrodinger particle with initial condition 
'f'o(C) ~ and 2 ^r-sH^—t) is its backward 

counterpart, with initial condition 4'o(C) + 

The evolution of the solution of the free particle KG 
equation is given in Fig. 21 where we have used an 
initially Macdonald distribution, namely for 4'o(fc) = 
^-p^/i+k'-^ ^i(^) = 0. Then, Eq. (IMl) gives 

'^KG{^,T) = [^j=b^Ki{c*) + ^^bK^(^c)\, (55) 

Lv27rc v27rc J 

where c and its complex conjugation c* are defined be¬ 
low Eq. (071). This is the effect of the interference with 
the backward time solution which yields a fairly rich evo¬ 
lution according to which the wave function undergoes a 
free propagation. The use of a Gaussian does not provide 
any significant difference. 

The evolution operator associated with Eq. (251) . 
related to two components Schrodinger-like equation. 


should be written as the matrix 2x2: 


Ukg{t) 


cos(t£)) sin(rL>)L> ^ \ 
— sin(rT>)T> cos(tZ)) / 


(56) 


whose determinant is equal to one. The inverse operator 
of Ukg{t) exists and it fulfils the relation \UKGij)\~^ = 
Ukg{—t)- In consequence Eq. (l48l) can be rewritten as 


/'I'o(e,r)\ 

Ui(e,r)y' 


Ukg 


(^om 

Ui(e)y ’ 


(57) 


where the equation on the upper component 
is given in Eq. (IMl) . whereas 4'i(^,r) = dr'^KGi^^T)- 
Moreover, the conserved quantity associated with such 
an evolution operator is /((^(l'I'o(C)P + I^~^l'I'i(C)P)'IC! 
which cannot be understood in terms of probability 
rather it has been constructed in analogy with the en¬ 
ergy of an harmonic oscillator. 

Let us observe that in Fourier space Eq. (IMl) reads 


bKG{k,T) = cos[Ta;(fc)]4'o(fc) + (58) 


where U!{k) = y'l + k'^ and 4'i(fc), i = 0,1 are the initial 
condition written in the Fourier space. The analogy with 
harmonic oscillator is made more clear by noting that in 
the Fourier space, the quantity 



[drbKGik,T)f + [Lu{k)bKGik,T)]'^dk 


(59) 


is a constant of motion. 

Finally, we point out that it is also possible express 
the solution of our problem in terms of the free particle 
Dirac solution by setting 

D = -iad^ + /3. (60) 


VI. QUANTUM RELATIVISTIC PARTICLES IN 
A LINEAR POTENTIAL 

A. The Salpeter equation with a linear potential 

Let us consider the Salpeter equation with a potential 
function U(^) 

idr'S{^,T) = y/f - (9|4'(C,t) -bMoI^(O) 

^(e,0) = 4 / 0 ( 5 ), (61) 

which in the Fourier space reads 

idrb{k, t) = \/l + fc2^(fc, r) -I- p,oVb{k, r), (62) 

where V = V{idk) and /rp is a dimensionless constant. 
We note that, in general case, the presence of an oper¬ 
ator V, which doesn’t commute with the kinetic term, 
hampers the possibility of a straightforward analytical 














solution for Eq. (1621) , which can however be obtained for 
V (^) being a linear function whose physically meaningful 
example is when a relativistic, charged particle interacts 
with a static constant electric field. In this case, Eq. dsa) 
can be rewritten in the form 

idr'^(k, t) = \/l + t) + i/ro9fc^(fc, t). (63) 

To solve Eq. (1631) first we eliminate the derivative by the 
use of the transformation 

^(/c,r) = (64) 


which yields a modified equation: 

i9.r4’(fc, r) = it'$(fc, r) (65) 

k{k, t) = y'l + pgTMoO. _ (00) 

The use of the procedure touched upon in Ref. [l^ allows 
to cast its formal solution in the form 


<I’(fc, t) = exp 


-i / k{k,T')dT' 
Jo 


^0 (fc) (67) 


where <I’o(fc) = ^o(fc)- The solution can finally be ob¬ 
tained in the form 


4'(fc, r) = r) = $(fc -|- /rpr, r) 

= e-*/oVl+(fe+^or'Pdr'^^(^ ^ ^^^)_ (0g) 

In Figs. [5] we report a comparison between the solu¬ 
tions for the free particle case and that of Eq. (1611) with a 
linear potential put in. In both cases we have assumed as 
initial wave function the function given in Eq. (1161) with 
/3 = 1. It is important to emphasize that the presence 
of the potential provides a suppression of the emergence 
of the bi-modal behaviour by inducing a kind of external 
localization. 


B. The Dirac equation with a linear potential 

The one dimensional Dirac equation with a linear po¬ 
tential has the form 


idr^ = —iad^^ + -I- (69) 


where fio = 


1 0 
0 1 


The solution of the problem is 


slightly more complicated than in the case of the Salpeter 
equation since the Hamiltonian is now expressed as linear 
combination of Pauli matrices. The formal solution of 
Eq. (l69l) can be given by using the operational method 
with the evolution operator dehned as 


_ ^Tad^+i^Tfto+iTP 


(70) 



FIG. 5. (a) Comparison between modulus of wave function 
evolution in presence of linear potential (I; blue curve; Eq. 
(1681) with ^{k) given in Eq. (1161) 1 and the wave function 
evolution of free particle (II; red curve; Eq. (1271) 1 at r = 4.3 
, (b) the same as (a) but r = 2.3 . 


commutator expansion and write the operator U in 
the following disentangled form 


u = g'i'iko-r+iT$^T&d^ ^i — fJ.oa 


(71) 


which is valid at the order The repeated applica¬ 

tion of the evolution operator dni on the initial state, 
yields the solution in the form of the following recursion 


3l„(e,T) 


giT(l-|-/ioC) 

0 


e 


— IT 


0 ^ 

(l-A^o?) J 


X 


cosh(r9j) sinh(r9{) 
sinh(rcl{) cosh(r9j) 


4L„_i(e,T), (72) 


where 






U-iOJo’ 


(73) 

(74) 


i’tAO = ^bP±{^ + T)±^±k-T)]. 


Albeit the problem is amenable for an exact treatment, 
we use the Zassenhaus formula at the second order in the 


with 


(75) 
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The previous iteration can be extended to any type of 
potential and specific relevant examples will be discussed 
in a paper specially dedicated to this problem. 


VII. CONCLUSION 

In this paper we have presented a number of issues also 
of practical nature to deal with problems concerning rela¬ 
tivistic equations and other fractional evolution problems 
often encountered in physics. We have left open many 
points and in particular the link of the present treatment 
with previous and well established methods. In particu¬ 
lar the use of Green function technique for the treatment 
of the Klein Gordon equation or the method associated 
with the use of Wightman solution . In this last case 
the solution of the free particle Klein Gordon writes, in 


terms of Wightman (and anti-Wightman as well) func¬ 
tions expressed as combinations of Hankel and Macdon¬ 
ald functions of the first kind. The relevant link with our 
formalism and with other type of solutions (Pauli-Jordan, 
Feynman, Stueckelberg) will be the topic of a forthcom¬ 
ing investigation. We have touched on many aspects of 
a relevant formalism and we have fixed the main math¬ 
ematical apparatus which can be exploited to discuss a 
number of problems in relativistic nonlocal quantum me¬ 
chanics. 
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